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PREFACE 



The Thermal Pollution Croup at the University of Miami has been 
developing three-dimensional mathematical models for predicting the 
hydrothermal behaviors of bodies of water subjected to a heated effluent. 
Generally speaking, thase models can be classified in two categories, 
namely, the free surface models which take into account the variation 
in water surface elevations and rigid-lid models which treat the water 
surface as flat. 

To enable the prospective user to use these models as accurate 
predictive tools, particularly to assess the environmental impact in the 
case of cooling lakes, they have to be calibrated and verified at a num- 
ber of sites. The present volume describes the application of the rigid- 
lid model developed by this group to a rather complicated site, namely. 
Lake Keowee in South Carolina. Lake Keowee is rather unique since 
it is used by a Nuclear Power Plant as a cooling lake as well as two 
other hydroelectric stations which use it as lower and upper ponds. 

This is the final verification of these models concluding a series of 

such verifications made possible by funding and technical assistance 
provided by the National Aeronautics and Space Administration (NASA) 
and the Environmental Protection Agency (EPA). 

This model will eventually be made available to all prospective 
users by NASA and EPA. The present volume together with the 
"Three Dimensional Rigid-Ud Model User's Manual" is intended for 
assisting such future users. 



ABSTRACT 



The Rigid Lid was developed by the University of Miami, Thermal 
Pollution Croup, to predict three-dimensional temperature and velocity 
distributions in lakes. This model was verified at various sites (Lake 
Belews, Biscayne Bay, etc.) and the verification at Lake Keowee was 
the last of these series of verification runs. 

The verification at Lake Keowee included the following phases of 
work. 

1. Selecting the domain of interest, grid systems and comparing 
the preliminary results with archival data. 

2. Obtaining actual ground truth and infrared scanner data both for 
summer and winter. 

3. Using the model to predict the measured data for the above periods 
and comparing the predicted results with the actual data. 

The model results have compared well with measured data. Thus, 
the.model can be used as an effective predictive tool for future sites. 
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SECTION 1 
INTRODUCTION 



BACKGROUND 

Understanding the environmental impact of hot water discharge from 
the condenser of a steam power plant is of considerable importance in 
helping to preserve aquatic life and also the overall efficiency of the 
power plant. In recent years, with the improvement of computing facili- 
ties, a number of different typ^s of mathematical models have been de- 
veloped to predict such effects. These models. If property calibrated, 
are extremely useful in determining the complete environmental impact 
arising d'-ie to normal as well as abnormal operating conditions, which 
may be difficult to realise experimentally. This is particularly important 
with respect to estimating these adverse effects during the planning 
stages, thus, helping to minimize them in the final construction. 

The mathematical models solve the basic fluid flow and energy 
transfer equations subject to certain simplifications and assumptions. 
These may be used to determine the long term heat budget for the cool- 
ing lake as well as determine the detailed velocity and temperature dis- 
tributions within it. The accuracy of predictions using these models is 
extremely important and the only way to verify this is to apply these models 
to certain known sites and compare the predicted values with actual 
measured quanJties. Measurements may be obtained by measuring the 
water velocities and temperatures with anemometers and thermometers 
as well as by infrared remote sensing photography. Each of these 
techniques has its own advantages and disadvantages. To get a proper 
representative data base, all of these have to be used simultaneously. 
Other factors which affect the water conditions are the meteorological 
conditions, viz. wind speed, solar radiation, air temperature and humidity. 

OBJECTIVES OF PRESENT WORK 

For the past several years the thermal pollution group at the Univer- 
sity of Miami has developed a number of mathematical models for both long 
and short term simulation. The model considered here is a three-dimen- 
sional rigid-lid model designed for short time prediction of detailed three- 
dimensional velocity and temperature profiles in th*. region of the thermal 
plume. This model has been applied '-i the past years to several sites, 
namely, Lake Belews (Lee, Sengupta, Mathavan, 1977). The predictions of the 
model at the above sites compared reasonably well with measured data. 
The results of these runs heve been pub t shed by the University of Miami in 



•artier volumes In cooperation with NASA. 

The success of this model at the above sites led to the belief that 
this model could be used in future predictions. It was mainly for this 
reason that a site was chosen which was quite complex in nature and 
which would bring out the prediction capabilities and limitations of the 
model. This would also be of great help to future users of the model 
in order to obtain a proper interpretation of the results when used as 
a predictive tool. 

Lake Keowee in South Carolina provides such a site; hence it was 
chosen for the final verification of the model. A description of the site 
is given in the following section. 

DESCRIPTION OF LAKE KEOWEE SITE 

Lake Keowee is located on the north of the state of South Carolina 
about 40 km west of Greenville. It was made from 1968 through 1971 by 
damming the Little River and Keowee River. At present, it constitutes 
Ouke Power Company's Keowee Toxaway Complex. 

The lake has two arms connected hy a canal (maximum depth 30. 5 m). 
There are three power plants on it, nam«ly, the Oconee Nuclear Station, 
Keowee Hydro Station and Jocassee-Pumped Storage Station. The Oconee 
Nuclear Station is a three-unit steam-etectric station with an installed 
capacity of generating 2580 MW. The Oconee Nuclear Station draws in 
condenser cooling water from the lower arm of Lake Keowee and discharges 
the heated effluent to the upper arm of the lake. The intake structure 
for the condenser cooling water allows water from 20 to 27 m depth (full 
pond) to pass through. The discharge structure Has openings from 9 to 
12 meters deep through which the CCW returns directly to the lower arm 
of Lake Keowee. 

Lake J ocas see is located north of Lake Keowee and is used as a 
reservoir for Jocassee-pumped storage station. The upper arm of Lake 
Keowee also serves as the tower pond for Jocassee-pumped storage sta- 
tion. The Jocassee Station has reversible turbines with a maximum gene- 
rating flow (into Lake Keowee) of about 820 m'/sec and a maximum pump- 
ing flow of about 775 m'/sec; the net flow into Lake Keowee from Lake 
Jocassee is about 15.5 m'/sec. 

Lake Keowee has a full pond elevation of 243. 8 m above mean sea 
level. At full pond it has a volume of approximately 1.18 x 10 m 1 , an 
area of 74 km 2 , a mean depth of 15. 8 m and a shoreline of about 480 km. 
The outflow from Lake Keowee is through Keowee hydro station. The 
flow through Keowee hydro station varies from 1.4 m I /sec (leakage) to 
560 m J /sec (during peak load operation). The maximum draw down of 
the lake is 7. 6 m. 

A map of Lake Keowee is shown in Figure 1. The above data was 
obtained from Duke Power Company 1976. 



SECTION 2 
CONCLUSIONS 



The main objective of this work was to prepare a mathematical 
model package to assess environmental impact due to a heated effluent. 
Hence, the f inal test for the model was at Lake Keowee. The reasonably 
good agreement of the predictions with data was shown in the last sec- 
tion. The siape of the plume, as predicted by the model resembles 
closely but does not correspond exactly to the plumes obtained from the 
IR scanner photographs. An important comparison in this case is the 
area covered by each isotherm which gives an indication of the spread 
of the heated effluent. This area is approximately the same in the mea- 
sured and predicted cases for both summer and winter runs. Another way 
the accuracy of the model was determined was by the root mean square 
deviations of the temperature between predicted and measured values 
over the entire domain. These deviations are shown in Table 1U. 

The accuracy of the IR scanner isotherms is about 0. 5 C C (from the 
sensitivity of the process). The accuracy of the position of the isotherms 
is within ±0. 5 x (grid space) in this case resulting from the lateral dis- 
tortion in the Digicolor map and the process of transfering this on to 
the computational grid for the purpose of comparison. 

Based on the above results, it may be concluded that the predictions 
mac' by the model are reasonable beyond doubt when applied to cooling 
lakes or a similar site. Hence, the model may be used as a predictive 
tool for obtaining three-dimensional temperature and velocity profiles in 
the vicinity of a thermal effluent discharge, and the results may be used 
in evaluating the performance of existing or future thermal power plants. 



SECTION 3 
RECOMMENDATIONS 



Various numerical models have been developed to study the effects 
of heated discharge and meteorological conditions on bodies of water. 
Most of these models are one or two dimensional. These models have a 
high computational speed but only give horizontally or vertically averaged 
values of temperatures. 

Three-dimensional models, however, have a much finer resolution 
but they consume larger computer time. The three-dimensional rigtd- 
lid model can be used to obtain detailed temperature and velocity distri- 
butions in a domain where surface gravity waves are small compared to 
the depth of the domain. This model, as compared to free-surface mo- 
dels, runs faster since surface gravity waves are eliminated by this 
rigid-lid assumption. 

A proper method of using this model would be to run a one-dimen- 
sional model initially to obtain a rough picture of the temperatures and 
then using this model to obtain a better resolution, the 1-D results being 
used as ambient conditions. 

The following improvements have been suggested for the model. 

1. Since all natural flows are turbulent, proper turbulent closures are 
needed to make the model meaningful. At present, the simplest 
possible closures, namely constant eddy viscosities ani c-ddy diffu- 
sivities, have been used. However, better results may be obtained 
by using a higher order closure. 

2. At present, the model uses uniform horizontal grids and stretched 
vertical grids. Nonuniform horizontal grids could be introduced for 
better resolution near the boundaries. 

3. The program has been written to be run as a batch- job on the com- 
puter. It could be made interactive so as to enable the user to run 
it on a terminal. However, this would require some modifications in 
order to reduce the storage space. 



SECTION « 
MATHEMATICAL FORMULATION AND MODEL DESCRIPTION 

CHOICE OF MODEL 

Two three-dimensional models have been developed by the thermal 
pollution group at the University of Miami. The first one is a free-sur- 
face model, which takes into account the variation in surface heights of 
a basin. The second one is the rigid-lid model which assumes the verti- 
cal velocities on the water surface as zero. These models have been de- 
scribed in previous publications of this group. The choice of a particular 
model for simulating flows depends on the nature of the site. 

Lake Keowee is a relatively small, closed basin with maximum depth 
of 30 m. The small horizontal numerical grid dimensions demand extremely 
small integration time steps to satisfy the Courant-Lewy-Fredrichs condi- 
tions inherent in free-surface formulations. Since the surface waves are 
small compared to the depth, i.e. h/H<<1, a rigid-lid formulation is suit- 
able. This rigid-lid model allows acceptable accuracy with acceptable 
time step size since the C.L.F. criterion is eliminated. This makes the 
rigid-lid model a rather obvious choice. 

The rigid-lid model has been applied previously to Biscayne Bay and 
Lake Belews sites with reasonable accuracy, as mentioned in the previous 
section. This led to the belief that it would be ideally suited for the 
Keowee site. 

DESCRIPTION OF THE MODEL 

(Portions of the following section have been published earlier by 
this group.) 

The rigid-lid model has the following capabilities: 

1. Predict the wind-driven circulation. 

2. Predict the circulation caused by inflows and outflows to the domain. 

3. Predict the thermal effects in the domain. 

4. Combine the aforementioned processes. 

The model solves equations for fluid flow (momentum and continuity) 



and heat transfer together with the equation of state in a three-dimen- 
sional domain. Since most geostrophic flows are turbulent, the resulting 
Reynolds stress terms In the governing equations are replaced by eddy 
transport coefficients. The fluid is considered to be an incompressible 
Boussinesq fluid and the hydrostatic assumption has been made. 

The rigid-lid assumption imposes a zero vertical velocity condition 
on the surface without affecting the horizontal velocities. This causes 
the surface pressure to be different from atmospheric which under special 
conditions may be translated into surface gravity wave heights if no lid 
is present. This assumption distorts transient time scales but does not 
affect circulation patterns as has been demonstrated by Berdahl (1970), 
Crowley (1969, 1970), Hag and Uck (1973) and Young, Liggett and 
Gallagher (1976). Since surface gravity waves no longer have to be 
reproduced, computer time is saved considerably and this assumption is 
quite adequate for cooling lake studies where surface waves are not 
large and do not change rapidly. 

To convert regions having irregular bottom topography to constant 
depth regions for computation stability a vertical stretching, suggested 
by Freeman et al. (1972), has been used. 

GOVERNING EQUATIONS 

The vertical stretching used in the rigid-lid model, originally 
adopted by Sengupta and Lick (1974), is of the form 



Z 



h(x,y) 

This enables the same number of grid points to be used all over 
the domain without using variable grid sizes. The resultinq nondimen- 
sion?l governing equations in the new (a, 8,y) coordinate system instead 
of ix.y.i) are 

Continuity Equation: 

li»2idi + ilUlI + h iiL, o 

3a 3S 3y 

Momentum Equation: 
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Hydrostatic Equation; 
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Energy Equation: 
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Equation of State: 

3 = 1.029431 - 0.000020 T - 0.0000048 T 2 

where 3 is in gm/cc and T is in °C 
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and w(2=0) = (rigid lid) 
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Quantities with subscript Yef are reference quantities; H and L are 
vertical and horizontal length scales. The variables with wavy lines on 
top are dimensional quantities. 
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If Prandtl number is equal to 1, then A , = B -. 

A u and A are the eddy viscosities in horizontal a 
n V 

B u and B are the eddy diffusivities in horizontal and vertical directions, 
n v 



A u and A are the eddy viscosities in horizontal and vertical directions 
rl V 



To obtain a predictive equation for pressure, the horizontal momentum 
equations are integrated from Z=0 to Z=h, where h is the nondimensional 
depth h/H. The integrated equations are then differentiated with respect 
to a and S summed. 

The Poisson equation for surface pressure becomes: 

3 2 P 3 2 P , 

3ot 3 3B* h 3a x, x. x p 

+ Z T-J-* "A +C-Y) 
h 3B l y. y, y p 



The last term is the Hirt and Harlow (1961) correction term which accounts 
for nonzero vertical velocities at the rigid lid. The variables (B , B , 
and A , A etc.) are given below: x ^ 

*1 X 2 
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The set of equations (2«*-30) together with appropriate boundary 
conditions constitute the mathematical model. The model has been de- 
scribed in earlier publications by this qroup (Lee and Sengupta, 1976). 

INITIAL AND BOUNDARY CONDITIONS 

The nature of the equations requires initial and boundary conditions 
to be specified. As the initial condition, the velocities, temperatures 
and densities are specified throughout the domain. Boundary conditions 
are specified at the air- water interface, geographical boundaries of the 
domain, the bottom of the basin and the efflux points. At the air-water 
interface a wind stress and a heat transfer coefficient are specified. 
The conditions on the lateral wails are no-slip and no-normal velocity 
for the momentum equations. These walls are assumed to be adiabatic. 
At the floor of the basin, the conditions of no-slip and no-normal velo- 
cities are applicable. The energy equation has a heat flux boundary 
condition. 

At points of efflux open-boundary conditions are specified. If the 
flows are known at the points of efflux, the flow velocities are specified. 
Usually, the temperature is known only at the point where the heated 
discharge (from the power plant) enters the domain. At points where 
the temperatures are not explicitly known, open-boundary condition (i.e. 
zero gradient condition) is specified. The same holds true if the velo- 
cities are not explicitly known (e.g. the connecting canal between the 



two arms of Lake Koewee) . 

Hence, the boundary conditions are: 
At the surface 

Y = 
fl = (Rigid Lid) 

3u f hH . 

3y " l U jT^zx 
' ret v 

, T hHK 

3 v , hH A . 

where t and t are wind stresses in the x and y-directions respectively. 

T_ is the equilibrium temperature. 

T is the water surface temperature. 

K is the surface heat transfer coefficient. 

t and t are calculated from the wind velocity usinq Wilson 
zx zy curve (B. W. Wilson, 1960). 

K and T are calculated as follows: 
s e 

K =U.5 + 0.05T + BF(w) + 0. 47(w) (wrrT^C' 1 ) 

where T is in °C. 

F(w) =9.2 + 0.46 W 1 (wm'^mHg" 1 ). 

B = 0.35 + 0.015 T + 0.0012 T *(mmHg/°C). 

T = (T * T All (°C). 
m so 

T . = Dewpoint temperature (°C). 

T c = T . + H /K . 
Ed s s 

H = Surface solar radiation (w/m 2 ). 
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At the bottom of th« basin 



On lateral walls 
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SPATIAL DIFFERENCE SCHEMES 

The numerical solution of the momentum and energy equations are 
explicit. The values of velocities and temperatures at a future time are 
determined completely using the values at the present and previous 
time steps. Finite difference forward time and central space difference 
schemes are used. Diffusion terms are written using a Dufort-Frankel 
finite difference scheme to relax the diffusive stability criterion. The 
convective stability criterion, hov/ever, is not affected. 

A horizontal staggering is used in the computational grids. Hori- 
zontal velocities and temperatures are calculated at the main- grid points 
while vertical velocities and pressures are calculated at half-grid points 

The predictive Poisson equations for calculating rigid-lid pressures 
is finite differenced using a five-point scheme. This is solved by 
successive over relaxation (Liebmann Method) . Terms on the right hand 
side of the pressure equation are obtained bv integrating terms in the 
horizontal momentum equations over the depth using the trapezoidal rule. 

At the boundaries, single-sided schemes are used since the boundary 
points do not have two adjacent points. A curve is fitted through the 
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two most adjacent points towards the interior of the domain. The values 
of the variables, where they are not defined, are obtained by averaging 
the values at four points around the point where the values are defined. 

STABILITY 

It is not possible to make a strict stability analysis of the system 
of equations under consideration. It is customary, however, to take 
advantage of the stability analysts for the one-dimensional Burger's 
equation since this contains an unsteady term, a convective term and a 
diffusion term. 

In the present case these criteria can be written as follows for the 
choice of the time step At 

Convective , r „ ii 

At < -rj- 

where U is the maximum horizontal velocity in the domain. 

Diffussive . f . (Ax) 2 

MARKER MATRICES 

Since natural bodies of water have irregular boundaries, the physical 
boundaries have to be approximated in a rectangular coordinate system 
using marker matrices (Fortran variable MAR for the full grids and MRH 
for the half grids). The convention used is as follows: 

MAR = 0, point outside the region of interest (i.e., on dry land) 

In approximating such a boundary using a rectangular grid system 
only a portion of the resulting grid falls within the water. To prevent 
computations to be carried out on dry land markers have to be used to 
distinguish between points lying within and outside the region of in- 
terests. Such markers have to be used both for the main and half-grid 
points. Fortran symbols used are MAR for the main grid and MRH for 
the half-grid system. The convention used is as follows: 

MAR = 0, point outside the region of interest (i.e., on dry land). 

MAR = 1, point on far y-boundary. 

MAR = 2, point on near y-boundary, 

MAR = 3, point on near x-boundary. 

MAR = 4, point on far x-boundary. 
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MAR - S, outside corner, on near x-boundary and far y-boundary. 
MAR = 6, inside corner on far x-boundary and far y-boundary. 
MAR = 7, outside corner on near x-boundary and near y-boundary. 
MAR = 8, inside corner on near x-boundary and near y-boundary. 
MAR = 9, outside corner on far x-boundary and near y-boundary. 
MAR - 10, outside corner on far x-boundary and far y-boundary. 
MAR =11, interior points (within region of interest). 

Similarly for the half-grid points. 
MRH = 1, corner on far x-boundary and far y-boundary. 
MRH = 2, points on near y-boundary. 
MRH = 3, points on near x-boundary. 
MRH = 4, corner at near x- and near y-boundaries. 
MRH - 6, far corner on x-axis. 
MRH = 7, corner at far x- and y-boundaries. 
MRH = 9, interior grid points. 
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SECTION 5 
APPLICATION TO LAKE KEOWEE 



INTRODUCTION 

The three-dimensional rigid-lid model described in Chapter II has 
been applied to Lake Keowee, South Carolina, to predict the wind-driven 
circulation, circulation caused by inflows and outflows, and the thermal 
dynamics of a region of this lake. This region includes the buoyant 
plume, the Keowee hydro dam, the connecting canal and an open-boun- 
dary where the effects of inflow and outflow to Jocassee-pumped storage 
station are felt. 

Lake Keowee is a warm monomictic lake having one circulation period 
during the year, beginning in fall, and is thermally stratified during the 
summer. During the circulation period, the. water mass is vertically 
mixed. This type of lake can also be termed holomictic (Hutchinson, 
19S7). Lake Keowee could also be classified as a subtropical lake as 
it exhibits thermal stratification, a period of total circulation preceded 
by the fail overturn, and surface temperatures above 4°C. The highest 
surface temperatures typically occurred in July and August and ranqed 
from approximately 27 to 29°C. During the periods from November tn rough 
February the lake was isothermal or nearly so. It can be concluded from 
the above that a series of meteorological events primarily govern Lake 
Keowee's thermal characteristic. However, as will be shown below, the 
near field region which makes up the region of interest in this present 
study, is very strongly affected by the buoyant plume, inflow and out- 
flow and the Jocassee-pumped storage station. 

CHOICE OF DOMAIN AND GRID SYSTEM 

The region of interest has already been mentioned in the last sec- 
tion. Here, an attempt will be made to justify this choice and then a 
description of the grid system will be given. 

The study and prediction of the thermal impact on the hydro-dyna- 
mics of Lake Keowee has been done in two main sections: 

1. Using a 1-D continuous temperature profile prediction model over 
the period 1971-1979, the area of interest includes the two main 
arms of Lake Keowee and the connecting canal. The application of 
this model has been published by the thermal pollution group of the 
University of Miami (Sengupta et al 1980). 
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2. Using the 3-D rigid-lid model, which is described In this report. 

The division has been necessary since the <-0 model predicts 
temperature profiles continuously over a period of time (several years 
if required) ; it can be used to predict the initial conditions required to 
run the 3-D model. The 3-0 model, on the other hand, has better 
resolution and is therefore used for predicting the near-field thermal 
impact on the lake. This near-field area of interest is shown in Figure 
5. In selecting this area some preliminary runs were made. These are 
Ascribed later. 

The region of interest is divided into a square grid as shown in 
Figure 5. The I and J axes (corresponding to u and v- velocity axes) 
are numbered as shown in the above figure. For the top 'open boundary', 
'J' increases from 7 to 20. The left side boundary shown 'I* increasing 
from 1 to 17; the right side boundary shows 'I' increasing from 1 to 16, 
and finally, the bottom boundary shows 'J' increasinq from 1 to 18. The 
orientation of these boundaries with the north-south direction is also 
shown in this figure. The squares are (152.4 x 152.4) square meters. 
Arrows indicate directions of flow into or out of domain. 

SUMMARY OF DATA 

For the purposes of calibrating and verifying the model, an archival 
data base was established and two remote sensing and field data collection 
missions (summer and winter) were undertaken. 

Archival Data Base 

The archival data base are summarized in Figure 6 and Tables 1 
through 4. These figures were taken from Duke Power Company (1976). 
Figure 6 shows the measured surface isotherms for September 10, 1975. 
This figure constitutes the data base on which the comparisons for the 
preliminary and archival runs are based. The runs are described below. 
Table 1 shows the monthly gross thermal capacity factors for Oconee 
Nuclear Station (1973-1977). Table 2 is a summary of Oconee Nuclear 
Power Plant for September 10, 1975. Table 3 shows the input data 
Used for the preliminary and archival runs. The data include the flow- 
rate, inlet temperature, discharge temperature and velocity, ambient 
temperature, depth at discharge, discharge width, air temperature, 
wind speeds and the verticil and horizontal eddy diffusivities. Table 4 
is a summary of the volune and area data for Lake Keowee. 

August 1978 Data 

To obtain an understanding of seasonal behavior of Lake Keowee, 
data was gathered both in summer and winter. This data was used to 
verify temperatures predicted by the model. Two extreme weather 
conditions were selected with the idea that if the model predictions were 
accurate for these conditions it would be so for other intermediate con- 
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dltions (viz., spring or fill). 

(n spite of the fact that the modal was used for predictions in the 
vicinity of the thermal plume, temperature and velocity readlnqs were 
taken for the entire lake. This was done in order to ensure complete- 
ness of the data and also to check the choice of the domain of interest 
within the lake. 

Summer data missions were carried out during the period of August 
24 and 25, 1978. The data collection team consisted of representatives 
from NASA, EPA, University of Miami and Duke Power Company. Three 
boats provided by Duke Power Company were used in the collection of 
ground truth data. The data collection stations »rt shown in Kigure 7. 
This included measurements of water temperatures at various depths 
(from the surface to about 30 m) using YSI type thermistors and velocity 
measurements using an Endeco Type 110 current meter. Tha thermistors 
were calibrated before each set of readings using a mercury thermometer. 
At each measuring station the boat was anchored and the thermistors 
and velocity meters were lowered by cables. The cables were marked 
thus indicating the depth below the water surface where readings were 
token. The current meter indicated both the magnitude and direction 
of the horizontal component of the water velocity. The main problem 
encountered here was that most of the velocities were very small and 
close to the threshold of the instrument. Hence, drogues were used to 
determine average surface velocities. 

While the ground truth data was being collected, overflights were 
made by NASA aircraft to obtain synoptic IR scanner data. The aircraft, 
a twin engined (Beechcraft) (NASA-6), was specially equipped for this 
purpose. "* Flights were made at altitudes of 1000, 2000 and 3000 ft. A 
diagram of the flight plan is shown in Figure 8, The scanner photo- 
graphs were taken with window settings between 78 and 84°F. To cor- 
rect the IR scanner photographs with the ground truth data for each 
set of flights, the water surface temperature was measured with the 
aircraft overhead for at least one station. This was to correct the scan- 
ner readings due to errors caused by water vapor attenuation. 

The scanner used was a daedalus series DS-1250 and remote sensing 
of 8-14 um radiation is achieved by a Hg: Cd: detector. The detector 
had a 0.015 inch square sensitive »r^». which was optimum for the reso- 
lution and temperature sensitivitv required. This detector was mounted 
in an end-looking, metal cored cewar which had sufficient liquid nitrogen 
coolant capacity for approximately six hours of operation before refilling. 
Thic system p rejected through the floor of the airplane, NASA 6, and 
had a ":an angle of 120° centered above the vertical. A horizontally- 
mounted telescope with its axis along the direction of flight of the air- 
plane was contained within the scanner. A mirror rotating at 3600 RPM 
and mounted at 45° to the telescope directed heat radiation from the 
ground into the system. A one-third revolution of the mirror covered a 
complete step perpendicular to the scanner axis. Optical resolution 
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obtained by this method was about 1.7 milliradians, so the ground areas 
detected become a function of flight altitude; the data accuracy is 0. 5°C. 

The video signal from the infrared detector was amplified and re- 
corded on magnetic tape in the aircraft. A method developed by Daedalus 
called Digicolor was used to convert this stored information directly into 
color coded strip imagery. This process limits the number of output 
colors to light (from white 'hottest' to black 'coldest') for any input con- 
dition. The six colors between white and black indicate the six calibrated 
levels of the set of interest. The scanner's thermal reference sources 
were present in flight to 66°F and 8*»°F respectively, and the settings 
were recorded on the same track with the detector video to insure accu- 
rate voltage relationships irrespective of all amplifier gain adjustments. 

The color bands in the final Digicctor map indicate zones of con- 
stant temperature within the accuracy of resolution of the scanner; hence, 
the line formed by the junction of two adjacent color bands indicate an 
isotherm. The actual temperature of the isotherm is obtained by adding 
a ground truth correction term to the temperature indicated in the map. 

Altogether, three IR flights were made during the summer data 
collection mission. Out of these, two runs were made on August 24. 
The first run was from 0853 hrs EST. On August 25, a single run was 
made from 0908 hrs to 0953 hrs EST. 

The Digicolor maps were transfered on to enlarged maps of the region 
of interest so as to obtain a map of surface isotherms which could be com- 
pared directly with the values predicted by the model. Out of the three 
runs made, only two, namely August 24 morning and August 25 morning, 
were used for comparisons since the resolution of the colors in the re- 
maining Digicolor map was very poor. 

The values of Oconee Nuclear Stations flows and temperatures every 
hour are obtained from continuous water quality monitoring stations of 
Duke Power Company. Flow through J ocas see- pumped storage station 
and Keowee Hydro Station as well as meteorological conditions (viz., air 
temperature, wind speed and direction, humidity and incident solar radi- 
ation] were also obtained from continuous monitoring stations. The flows 
through Jocasse, Keowee and Oconee as well as the discharge temperatures 
are shown in Tables 5 and 6. The meteorological data (obtained hourly) 
are shown in Tables 7 and 8. Figures 9 and 10 show the hourly variation 
of Keowee Hydro Station flowrate on these two days. Figures 11 and 12 
show the same for Jocas see- pumped storage station. 

February 1979 Data 

The winter data mission was carried out during February 27, 28, 
1979. The same equipment were used in this mission. Three boats 
were used to measure water temperatures and velocities up to a depth of 
30 m. Two of the boats were equipped with thermistors only and were 
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used to measure temperatures in both branches of Lake Keowee. The 
third boat carried the current meter as well as a thermistor. The stations 
where readings were taken are indicated in Figure 13. Station #13 was 
used for ground truth correction of the IR scanner data. This point was 
chosen since it is virtually unaffected by the discharge, and the tempera- 
tures here remained fairly constant. 

Ground truth measurements were taken from 0950 hrs to 1353 hrs 
EST and 1600 hrs to 1833 hrs EST on February 27. Due to technical 
problems with the aircraft (NASA-6) the IR flights were delayed and 
were only from 1549 hrs to 1711 hrs EST on this date. The flight plan 
was identical to that used in the August data, mission. Black body set- 
tings of 38°F and 74°F were used and the flight altitudes were 2000, 
3000 and 1000 feet. 

On Feburary 28, 1979 ground truth measurements using three boats 
were taken from 0851 hrs to 1201 hrs EST and 1420 hrs to 1830 hrs EST. 
IR flights were run from 0850 hrs to 1002 hrs at altitudes of 2000 and 
3000 ft. The black body settings used were established by NASA. IR 
isotherms were constructed for the domain of interest at the University 
of Miami. These maps were used for verifying the results predicted by 
the computer. 

For obtaining meteorological data and flows through Oconee Nuclear 
Station, J ocas see- pumped storage station and Keowee Hydro Station, 
continuous monitoring stations of Duke Power Company, were used. 
Since the velocities in the lake were found to be extremely small, drogues 
were also used. The flows through the three power stations are shown 
in Tables 9 and 10. The data obtained was hourly and the variation of 
the flows through J ocas see and Keowee are shown in Figures 14 and 15 
and Figures 16 and 17 respectively. The meteorological data collected on 
February 27 and 28 are shown in Tables 11 and 12 respectively. 

The ground truth data collected during the summer and winter 
missions are presented by Duke Power Company (1978 and 1979). 

CALCULATION OF INPUT 

The 3-D rigid-lid model as described in Section 4 solves the three- 
dimensional momentum, continuity and energy equations. As shown be- 
fore, these equations are a set of nonlinear, coupled partial differential 
equations and require the following for complete solution. 

1. Initial values of the velocities and temperatures must be specified 
at all points within the domain. 

2. Boundary conditions for the above variables must be specified at 
all boundaries. 

The choice of the domain of interest as well as the grid system has 
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been described in earlier sections. At all solid boundaries the velocities 
and temperature gradients are specified as zero. At the surface the 
rigid-lid constraint of zero vertical velocities is assumed. The initial 
conditions assumed for starting the runs are zero velocities and con- 
stant temperatures everywhere within the domain. For subsequent runs 
the results of the previous run is used as the initial condition. Hence, 
the quantities yet to be specified are: 

1. Temperatures and flow velocities at the Oconee Nuclear Station dis- 
charge. 

2. Temperatures and velocities at the Keowee Hydro Station. 

3. The same at the Jocassee boundary and at the canal. (Items 1 and 2). 

4. Surface horizontal velocities and temperatures. 

The above quantities are termed as inputs to the model and the 
procedure used for obtaining them are discussed briefly below for the 
sake of completeness. For further details regarding the actual running 
of the programs, the reader is advised to refer to the 3-D riqid-lid 
User's Manual (Sengupta et al, 1980) prepared for this purpose. 
The calculations shown below were for February 27, 1979 simulations. 

Reference Quantities Used 

Reference length = L = maximum length of the domain = 2895.6 m. 

4/3 
Reference horizontal eddy viscosity A , = 0.002 L = 38311.48 cm 2 /sec. 

(Note: The constant •0.002' changes with different sites. This particu- 
lar value was used in running the model at Lake Belews and at Biscayne 
Bay. In this case the best value of the constant was found to be '.003* 
which yielded a value of A , = 60,000 cm 2 /sec. 

Reference depth = H = max depth considered = 16 m 

4 ; 3 
Reference vertical eddy viscosity A = 0.002 x H ' =37.43 cm 2 /sec. 

Reference velocity U . = 30 cm /sec. 

Reference time T - = LA/ , = 9652 sec. 
rer rer 

Oconee Nuclear Station Discharge Velocity 

The discharge is considered to take place through a point at a depth 
of 12 m (k=3). The discharge velocity is calculated as follows. 
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12 mi 




152. 1 m 



The total discharge = (100 ~ x V(|^) x 152.4 x 12) = Q 

(where Q = average discharge in m J /sec.) 

. \V = 7. 4207 cm /sec 

The average value of 'Q 1 over 24 hrs is taken since the variation 
is negligible. 



Nondimensional discharge velocity = 



Keowee Hydro Discharge Velocity 



ref 
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= 0.24740. 



The outflow through the Keowee hydro station is throuah a channel 
(152.4 m) x (12 m). 

The volume flowrate Q = (152.1 x 12 x V) m J /sec. 

(where V = discharge velocity (m/sec.) 

V = [Q/(152.4 x 12)] m/sec = i 52 . 4x^2x100 cm /sec 

Q is specified as a function of time with the help of polynomials and 
other functions. The curve is approximated and specified in subrou- 
tine INLET1; (refer to the user's manual (Sengupta et at, 1980). 

Keowee hydro flow approximation: refer to the user's manual (Senaupta et at , 1 980) 

(February 27 Data) 

SX = Conversion factor for converting discharge (incfs) to nondimensional 
velocity. 

SV1 = Nondimensional velocity. 

The velocities are approximated as follows: 

SV1 = 0.048 < TSDT < 6 
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SV1 ~SX * (((17.54- 0.048)72.) * (TSDT - 6.0) + 0.48) 

6 < TSDT < 8 
SV1 = SC * (((.048- 17.54J/4.) * (TSDT - 8.0) + 17.54) 

8 < TSDT < 12 

SV1 = SX * (0.048) 12 < TSDT < 24 

Jocassee Flow Velocity 

The entire flow to or from the J ocas see- pumped storage station is 
assumed to take place through the entire upper boundary. The flow 
through this area is shown in the following figure. 




The flow is assumed to be uniform over this area (i.e., equal flow 
velocities at all internal grid points within this area) and is assumed to 
take place simultaneously with the flow through Jocassee- pumped storage 
station. 

V = Q/K1S x 13 x 152.4) x 100] cm /sec 

(Q = flow through Jocassee-pumped storage station (m 3 /sec).) 

Q is positive when Jocassee is generating (i.e., the flow is into the 
region of interest) and negative when pumping (i.e., flow out of the 
region of interest) . 

Jocassee flow approximation: 

(February 27 Data) 

TSDT = Time from start of run (hrs) . 

SV = Velocity of flow through Jocassee boundary (nondimensional) . 

SF = Conversion factor to convert flowrate (cfs) to nondimensional 
velocity 

= 0.00322579 

The velocity is approximated as follows. 
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SV = SF * (-14.395 - (18.75 - 14.395) * (TSDT) 

1 < TSDT < 
SV a SF * (-18.754) 1 < TSDT < 5 

SV = SF * (((16.823 + 18. 7-54)7 3.). * (TSDT - 5.0) + 18.754) 

5 < TSDT < 8 
SV = SF * (((16.823- 0.1) /3.) * (TSDT - 8.0) - 16.823) 

8 < TSDT < 11 
SV =: SF * 0.1 11 < TSDT < 23 

SV = -SF * ((4.5 + 0.1) * (TSDT - 0.1) 

23 < TSDT < 24 

A similar procedure was followed for simulations of the other days. 

Condition at Open Boundaries 

Open boundaries are those where the values of temperatures and /or 
velocities cannot be specifically obtained but continuity of flow has to be 
maintained. One such boundary is at the mouth of the canal connecting 
the two arms of Lake Keowee. The condition — = is sepcified here 
both for velocities and temperatures. At the Jocassee and Keowee boun- 
daries, the same kind of zero gradient conditions is specified for tempera- 
tures only. The calculation of parameters (T^, TAU, etc.) for the sur- 
face boundary conditions are shown in Section 2. 



22 



SECTION 6 
RESULTS AND DISCUSSIONS 



The results of simulation using the 3-D rigid model at Lake Belews 
and Biscayne Bay are shown in previous publications by the University 
of Miami (Lee, Sengupta and Mathavan, 1977, end Lee and Sengupta, 
1976). In the above two cases the predictions made by the model agreed 
closely with IR scanner and ground truth data. The following section 
discusses the runs made for Lake Keowee. 

Preliminary Runs and Results 

These runs were primarily made for the selection of the boundary 
conditions, initial conditions and to test the behavior of the model when 
incomplete and /or arbitrary data are used. These runs are summarized 
in Table 13. 

1. Run Number L001 

(The Number L001 is a sequence number and must only be interpreted 
»s such.) 

This run was essentially made for debugging the computer program 
modifications. The features include a discharge velocity of 5.65 cm/ 
sec, a discharge temperature of 32. 3°C and a 16 meter constant depth 
region of interest. The total simulated time was 20.6 hours. The 
results are summarized in Figures 18, 19 and 20. The surface 
velocities after 8.64 hours are shown in Figure 18, while Figure 19 
shows the velocities in a vertical plant (1=11) after 21.6 hours. 
Since the effects of wind were not included in this run and there 
were no effects of J ocas see, the velocities are by the discharge 
velocity and the buoyant plume. The isotherm comparison of the 
archival run and this simulation run is shown in Figure 20. The 
ambient temperature is 29°C. The archival isotherms are higher 
than the predicted isotherms. This is expected since this run 
(without wind) was not undertaken primarily for comparison pur- 
poses, and the simulation did not reach st-.-Jdy state. 

2. Run Number L002 

This run is similar to 'LOOT but with variable depth of the region 
of interest. The simulation time was only 8.64 hrs. This time is 
long enough to determine how well the model handles a variable 
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depth domain. The surface velocities and vertical velocities (at 
1=11) are shown respectively in Figures 21 and 22. For complete- 
ness, the simulated surface isotherms are compared with the archival 
isotherms. This comparison is shown in Figures 23. For the same 
reason described in the previous subsection, the comparison cannot 
be expected to be more precise. 

3. Run Number LOO 3 

Run L002 is now repeated with the effects of wind included. The 
total time simulated was only 8.00 hrs. The reason for the short 
simulation time is as explained above . The surface and 
vertical velocities (J =7) are shown in Figures 24 and 25 respectively. 
A comparison of Figures 21 and 24 shows the small effects of wind on 
the surface velocities after such a small simulation time. The wind 
speed is 4 m/sec. A comparison of the predicted and archival iso- 
therms is shown in Figure 26. Compared with the other isotherm 
comparisons, this figure appears to have improved the difference be- 
tween measured and predicted. 

Archival Runs and Results 

Two runs were made (L004 and L005) using the archival data base 
described in the last section. 

1. Boundary Conditions 

The following specifications were used at the main input boundaries. 

a. At 1=1, J=8 to 19 for all depths (K) 

The water velocity at the Jocassee end is specified as 6. 9 cm/ 
sec. Open boundary is specified for temperatures as this end. 

b. At 1=11. J=1, K=3 

This discharge velocity is specified as 5.65 cm /sec, and the 
dishcarge temperature also specified as 32. 4°C. 

c. Finally, at 1=13, J=7, K=1 to 3 

Keowee dam discharge velocity is specified as 10.4 cm/sec. 

2. Run Number L004 

This r> the first archival run and was made with variable depth 
topog aphy for the near-field region of interest. 

3. Run Number L005 
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Run L004 is repeated but for constant depth region of interest. 
The domain was cut off at 16 meters from the water surface. This 
is the depth of the thermociine for Lake Keowee. 

4. Archival Results 

Figure 27 shows the comparison between the measured (archival 
; 9/10/75) and calculated (L004) isotherms. The ambient temperature 
was 29°C. This comparison can be seen to have improved when 
compared to those of the preliminary runs already described. The 
surface velocities are simitar *.o 'L005' which are shown in Figure 28. 
Flow through the Keowee drm can be seen very clearly on the right 
boundary of the figure, 'figure 29 shows the velocity field at J =7, 
which also shows the flow through the Keowee dam. The isotherm 
comparisons of measured (9/10/75) and predicted (L005) is shown 
in Figure 30. This is the best comparison of the runs described so 
far. The agreement is particularly good for the 29. 5 and 30°C iso- 
therms. The shapes of the measured and predicted 30. 5°C isotherms 
are similar. However, the measured area under the isotherm appears 
bigger. This discrepancy can be explained in terms of the nearness 
of this isotherm to the discharge. The total time simulated was 32.4 
hrs, which is also the time when steady state was reached. The 
temperature profiles for locations 1=11, J=7, K=1 to 4 through the 
simulated period are shown in Figure 31. The temperature profiles 
are vertical as expected since computations were carried out above 
the thermociine. A comparison of the surface velocities at this loca- 
tion is given in Figures 32. It can be seen from this figure that the 
u-component velocities predominate. It can also be seen from this 
figure that all the velocities stabilize after about 30 hours, showing 
steady state conditions. A similar plot. Figure 33, for the same 
variables but for location 1=11, J=2 and K=i is shown. Similar con- 
clusions hold for this figure. 

L001 through L005 were archival runs. The following section dis- 
cusses the results of simulations for the August data base (L006) 
and February data base (L007). Both these runs were made ! n 
stabs of 24 hrs of simulation at a time for a total period of 48 hours 
of simulation in each case. Results were stored and printed at the 
end of every hour of real time. The results at the end of the first 
24 hrs simulation were used as initial conditions for the next 24 hrs. 

Summer Runs 

Simulation was started from 0000 hrs August 24. The model was 
started 'cold* i.e., initial conditions used were zero velocities and con- 
stant temperature (29°C) throughout the domain. This conditions re- 
quired that the model had to be run for sometime before the effects of 
the cold start could be ignored. 

The inputs used in running the model have been discussed in the 
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previous section. The wind speed driving the upper layer varied from 
0.76 to 3.09 m/sec. The current or the flow through the Jocassee 
boundary varied from 0.61 to 1.12 cm /sec. The Oconee nuclear station 
discharge velocities and temperature were 7.42 cm /sec and 31.7°C re- 
spectively. (Average values were taken since the fluctuation was negli- 
gible.) All the above data were fed in every hour in running the model. 

The surface isotherms and velocities at each horizontal section 
(K=l, 2, 3, 4 and 5) were plotted using a Calcomp plotter. Results at 
the end of each hour were stored and plotted. The results used for 
comparison were the surface isotherms plots. Measured surface isotherms 
were obtained from IR scanner digicolor maps. By drawing both the 
data and the predicted values on the same grid system direct comparison 
could be made. 

The first IR runs were made on August 24 from 0853 hrs EST to 
1002 hrs EST. The isotherms drawn for this data are shown in Figure 
34. The isotherms shown are for 30. 5°C, 30.0°C and 29. 5°C. Only the 
above three isotherms fell within the domain of interest and represent 
the water surface temperatures at 1002 hrs EST, August 24, 1978. The 
results produced by the model are shown in Figure 35 and 36. Figure 35 
shows the same isotherm i.e., 30. 5°C, 30.0°C and 29. 5°C corresponding 
to 1007 hrs real time. It is seen that the temperatures predicted by the 
model are lower than the actual temperatures, showing that the plume 
spread was underpredicted. Figure 36 shows different predicted iso- 
therms corresponding to the same time. The temperatures of the iso- 
therms are 29. S°C, 29. 3°C and 29. 1°C. Comparing Figure 36 with 34 
it is seen that 29. 5°C in the IR data agrees very well with the 29. 1°C 
isotherm predicted by the model. The "measured 30. 5°C and 30.0°C iso- 
therms compared with the predicted 29. 5°C and 29. 3°C isotherms respec- 
tively. Hence, the errors in prediction in the three isotherms (30. 5°C, 
30. 0°C and 29. 5°C measured) are approximately 1°C, 0. 7°C and 0. 4°C 
respectively. 

The IR digicolor map for the August 24 afternoon data could not be 
used to construct isotherms because of bad resolution between different 
colors. Hence, the next useful comparison was made for the August 25 
morning data (0903-9553 hrs EST). The IR isotherms are shown in 
Figure 37. The same three temperatures, namely 30.S°C, 30. 0°C and 
29. 5°C. are represented here since the other isotherms lie beyond the 
domain of interest. Figure 38 shows the same isotherms as predicted by 
the model. The spread of the 29. 5°C isotherm is overpredicted while 
the 30.0°C isotherm is underpredicted. However, Figure 39 shows iso- 
therms for 29. 90°C, 29.70°C and 29,6°C as predicted by the model. 
These isotherms compare very well with the IR isotherms in Figure 37. 
The errors in this case are . 6°C for the 30. 5°C isotherm, 0. 3°C for 
the 30.0°C isotherm and -.1°C for the 29. 5°C isotherm. This shows a 
considerable improvement in the predictions as compared to August 2H, 
showing the effects of cold start gradually vanishing. 
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Figurt 10 shows the plot of the horizontal velocities of the surface 
as predicted by the model for 1007 hrs, August 24. During this time 
(1007 hrs EST), as can be seen from the Jocassee-pumped storage sta- 
tion flow and Keowee hydro station flow data, the flows through these 
stations were negligible. This leads to zero velocities at Jocassee boun- 
dary and zero velocity through the Keowee hydro discharge point in 
Figure 10. 

Figure 11 shows the surface horizontal velocities predicted by the 
model for August 25 (33.2 hrs run time). During this time (approxi- 
mately 1000 hrs, August 25) Jocassee-pumped storage station just started 
generating and Keowee hydro was generating. This accounts for the flows 
through the two boundaries. 

The main driving forces responsible for determining the shape of 
the isotherms are the ambient temperature, discharge temperature and 
flows through the Jocassee boundary. The wind is seen to affect the 
velocities in the upper layer only. This characteristic is displayed 
both in the data and the simulation results. Figure 31 (August 21 IR 
data) and Figure V (August 24 IR data). Looking at the Jocassee and 
Keowee flows it is seen that the flow through Jocassee is negligible dur- 
ing both these periods. The flow through Keowee is negligible at 10.00 
a.m. on August 24 but is 9352 cfs on 10 a.m., August 25. The average 
discharge temperature is also lower on August 25 as compared to August 
24. This causes the area under the isotherms in Figure 37 to be lower 
than those in Figure 34. The same difference is seen between the cor- 
responding predicted isotherms, namely Figure 36 (August 24, 1978) and 
Figure 39 (August 25) . 

Winter Runs 

Simulation was started from (0000 hrs) February 24, 1979. The 
model was started using zero velocities and constant temperature (= 10°C) 
as initial conditions. 

The inputs used for this run (L007) have been discussed in detail 
in previous sections. The Oconee discharge velocity (average value) was 
6.84 cm /sec and the discharge temperature (average value) was 18. 1°C. 
The wind speed varied from 1.61 to 4.52 m /sec and the Jocassee-pumped 
storage station boundary flow velocity ranged from 0.61 to 4.14 cm /sec. 

Values of velocities and temperatures were printed at the end of 
every hour. Surface isotherm plots and velocity plots were generated 
from these results. Figure 42 shows the IR scanner isotherms for 1648- 
1651 hrs, February 27. The temperatures are 13.0°C, 1 2. 5°C , 12.0°C, 
11.5°C and 11.0°C. The same isotherms us predicted by the model are 
shown in Figure 43. Ther* correspond to 17.12 hrs after the cold 
start. Comparing these two figures it is seen that the isotherms pre- 
dicted by the model have a greater spread. Figure 44 shows predicted 
isotherms corresponding to 13. 75°C, 13.0°C. 12.75°C, 1 2. 5°C and 12. 0°C. 
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This figure agrees better with Figure 12 showing an error of prediction 
of 0. 75°C, 0. 5°C, 0. 75°C and 1°C respectively for each of the isotherms. 

Figure 45 shows IR scanner isotherms for February 28, 1979 corre- 
sponding to 0948-0957 hrs E5T. The two isotherms within the domain 
are for 13. 0°C and 12. 5°C. Figure 46 shows the same two isotherms as 
predicted by the model. These correspond to 34.2 hrs of run time after 
the cold start at approximately 10 a.m., February 28. Figure 46 is in 
excellent agreement with Figure 45. Isotherms for IR data corresponding 
to February 28 afternoon could not ba drawn because of bad resolution 
among different colors in the digicolor map. 

Figure 47 shows velocities at the surface as predicted by the model 
for February 27, 17.12 hrs after the cold start. The flows through the 
Jocassee and Keowee hydro boundaries agree with the data obtained. 
Figure 48 shows the same for February 28, 34.2 hrs after cold start. 
Excellent agreement is obtained in this case between the velocities shown 
on the map and measured flows at Keowee hydro station and Jocassee- 
pumped storage station. As in one of the summer runs, the isotherms 
are affected by the Oconee discharge, Keowee hydro and Jocassee-pumped 
storage station. Figures 42 and 45 (isotherms for February 27 and 28 
respectively) show that the isotherms have spread out further on the 
second date. Keowee hydro and Jocassee-pumped storage station flows 
are negligible on both the dates. The average Oconee discharge tempera- 
ture during 1600 hrs, February 27, was lower than that during 1000 hrs, 
February 28. This accounts for the greater spread of the isotherms. 
The same is depicted in the results produced by the model. 
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Figure 4. MRH marker matrix (half grid points) 



. 1 1 U.UiU 



Di schatge 




\ 



Out. £ 1 *)w 
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Figure 6. Measured isotherms (Archival 9/10/7E) 
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Figure 8. Flight plans for IR data (Auoust and February missions) 
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Figure 9. Keowee hydro discharge data {Auqust 21, 1978) 
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Figure 10. Keowee hydro discharge data (August 25, 1978) 
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Fiaure 11. Jocassee-pumped storaae station discharae data 
(Auaust 24, 1978) 
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Figure 12. J ocas see- pumped storaae station data 
(August 25, 1978} 
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Fiqure 13. Keowee February 1979 mission showinq stations 
(CT = qroundtruth, FL3 = fliqht NoFL3) 
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Fiqure 14. J ocas see- pumped storaqe station discharae data 
(February 27, 1979} 
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Figure 15. J ocas see- pumped storaae station discharge data 
(February 28, 1979) 
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Fiqure 16. Keowee hydro discharae (February 27, 1979) 
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Figure 17. Keowee hydro discharge data (February 28, 1979) 
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Figure 18. Velocities at K = 1, Lake Keowee (rigid-lid model) 



47 



HUH K9» 1831 . 

DISCHARGE YEL3CITY i S.SSCM/SEC 

0I3CKR33E TEMPERATURE t 32.3'c 

HIN2 J O.QCH/SEC 

CURREHT » Q.O Cn/SEC 

T8TRL SlflULRTEO TIflE 1 21.6CHSS 



SCALES CH3SII3NTnLl 



0.00 
LENGTH SCfiLECflETERSJ I 



MlO 1 



51 .00 
i 



VEL3CITY SCALECCfl/SEC) °il£5. 



12. CO 



5CALE3 (VERTICAL) 



LENGTH SCflLEiMETEnSJ 



0.00 

i 



20.00 



YEL3CITY SCfiLECCfl/SEC) 



0.00 



5.0C 
i 




'S/1 



-* "V N * 



I \ \ » t I ^ • • » 



Figure 19. Velocities at I = 11, Lake Keowee (rigid-lid model) 



48 



^ 


RUN NQ : LQQ 1 • 

DISCHARGE VELOCITY 

DISCHARGE TEMPERATURE 

WIND 

CURRENT 

TOTAL SIMULATED TIME 

LENGTH SCALEC METERS) 

VELOCITY SCALECCI1/SECJ 


i 

5 . 6 5 C M . ' S E C 
32.3 a C 
O.QOM/StC 
Q.Q CM/SEC 
9 .S-l HRS 

*1Q' 
0-00 


61 .00 

1 


0.00 

1 


• 

30. GC 
j 



ARCHIVAL DATA i'o/i(}/7*5^ 
PREDICTED 



— .< 



s: 
























0..lV 


. ^ 


., 















r n.~ > 



•^ 



1 ^ 



Fiqure 20. Isotherms: measured and predicted (temperatures above ambient) 
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Figure 21. Velocities at K = 1, Lake Keowee (riqid-lid model) 
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Fiaure 22. Velocities at I = 11, Lake Keowee (riqid-lid model) 
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Figure 23. Isotherms: measured and predicted (temperatures above ambient) 
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Figure 2*». Velocities at K = 1, Lake Keowee (rigid-lid model) 
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Figure 25. Velocities at J = 7, Lake Keowee (rigid- lid model) 
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Figure 26. Isotherms: measured and predicted (temperatures above ambient) 
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Fiqure 27. Isotherms: measured and predicted (temperature* above ambient) 
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Figure 28. Velocities it K ■ t, Lake Kmwm (rigid-lid model) 
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Figure 29. Velocities at J « 7, Ukt Koowm ( rigid- 1 W model) 
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Figure 30. Isotherms: measured and predicted (temperatures above ambient) 
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Figure 31. LOOS, temparatura profiles 
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Figure 32. LOOS, velocity vs. time, I a 11, J ■ 7, K ■ 1 
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Figure 33. L005, velocity vs. time, I = 11, J = 2, K = 1 
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Figure 34. IR data cort««pondinq tQ 1002 hr«, Auqust 24, 1978, Lake Kaowaa 
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Figure 35. Isotherms at K = 1, Lake Kaowee (rigid-lid modal), 
simulations for August 24, 1978 
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Flqure 36. Isotherms at K ■ 1, Lake eowee (rigid-lid model), 
simulations for Auqust 29, 1978 
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Fiqure 37. IK data correspondinn to 0903-0953 hrs, Auaust 25, 1978 Lake Keowee 



66 





a 


RUN N 

Disca 

DISCH- 
WIND ! 
CURRE 
TOTAL 

LENGT 


0: LOOS 

fcRGE VELOCITY 
fcRGE TEMPERATURE 
3FEED (MAX) 
INTCJOCASSEE FLOW): 
. SIMULATED TIMS 

H SCALE (METERS) 


7.42 CM /SEC 
: 31.7 °C 
2.50 M/SEC 
2. 2 CM /SEC 
33.23 HRS 

10 1 
0.00 1 


91. ( 


)0 








T~ 


* 


















* 














































29.5 ^^^— ^^ 


■^^ 








.. 




\ 












) 


> 
1 
1 
































/ 


y 


" 
















30.0 






I 




fc » 







Figure 38. Isotherms at K « 1. Lake Keowee (rigid-lid model), 
simulations for Auaust 25, 1978 



67 



RUN NO: LOOS 
DISCHARGE VELOCITY 
DISCHARGE TEMPERATURE 
WIND SPEED (MAX) 
CURRENT (JOCASSBE FLOW) 
TOTAL SD4ULATED TIME 



7.42 CM /SEC 
31.7 °C 
2. SO M/SBC 
1.1 CM /SEC 
33. 2 HRS 



^ 



10 J 



LENGTH SCALE (METERS) 



0.00 

I— 



61.00 




Figure 39. Isotherms et K » 1. Leke Keowee (rigid-lid model), 
simulations for August 25, 1978 
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Figurt 40. Velocities it K * 1, Lake Ksowm (riqid-lid model), 
simulations for August 24, 1978 
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Figure «l. Velocities it K ■ 1, Lake K*owm (riqid-l!d model) , 
simulations for August 25, 1978 



70 




Figure 42. IR data corresponding to 1649-1651 hrs, February 27, 1979, 
Lake Keowee 
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FIgur* 43. \§oVmm* at K * 1, Lak« Ktow«« (Hpld-Hd mod«0, 
simulations for Fsbrutry 27, 1979 
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Figure 44. isotherms at K -1. Lake Keowee (rfqid-Hd model) 
simulations for February 27, 1979 
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Figura 15. IR data corresponding to 0948-0957 hrs, February 28, 1979 
Lake Keow^e 
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Figure 46. Isotherms at K =1, Lake Keowee (rigid-lid model), 
simulations for February 28, 1979 
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Figure 47. Velocities at K = 1, Lake Keowee (riqid-lld model), 
simulations for February 27, 1979 
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Figure 48. Velocities at K = 1, Lake Keowee (riqid-lide model), 
simulations for February 28, 1979 
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Table 1. Monthly Gross Thermal Capacity Factors (Percent) 
for Oconee Nuclear Station . 



MONTH 


1973 


1974 


1975 


1976 


1977 


Jan. 





13 


37 


30 


90 


Feb. 





28 


22 


49 


77 


Mar. 





29 


57 


47 


72 


Apr. 





26 


61 


27 


92 


May 


4 


5 


76 


33 


78 


Jun. 


9 


43 


78 


64 


46 


Jul. 


15 


38 


80 


59 


44 


Aug. 


15 


31 


82 


90 


36 


Sep. 


16 


46 


76 


77 




Oct. 


23 


19 


75 


61 




Nov. 


23 


31 


92 


52 




Dec. 


27 


30 


92 


55 




Annual 
Average 


11 


28 


69 


59 





Data based on three* unit capability for entire period. 

nllVL^^tr Actual MWH(t) X 100 

for ONS(%) " 7704 MW(t) x Hours in Period 
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Table 2. Operating Conditions of Oconee Nuclear Power Plant 
at Lake Keowee (September 10, 1975) 



Plow Rate 

Inlet Temperature 

Outlet Temperature 



7087 m/rain 
28.8 °C 
32. 4 °C 



The hourly averages of wind, air temperature and 
relative humidity: 



Time 


Wind 
Speed (MPH) 


Wind 
Direction (°N) 


Air 
Temp. 


Relative 
Humidity 


9:00 


7.7 


15 


69.0 


100 


10:00 


6.8 


65 


74.0 


70 


11:00 


8.3 


80 


77.0 


60 


12:00 


7.2 


95 


77.0 


63 


13:00 


9.7 


115 


78.0 


88 


14:00 


11.3 


205 


79.0 


84 


15:00 


- 


- 


77.0 


66 


16:00 


10.2 


75 


73.0 


84 
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Table 3. Input Data for the 3-D Model (Lake Keowee) 



Flow Rate: 


118.12 m'/sec 


Inlet Temperature: 


25.6 °C 


Discharge Temperature: 


32. 4 °C 


Discharge Velocity: 


5. 65 cm /sec 


Ambient Temperature: 


28. 9 °C 


Depth at Discharge: 


30 meters 


Discharge Width: 


152. 4 meters 


Air Temperature: 


28.11 °C 


Wind Speed: 


12.83 ft /sec 


Vertical Eddy Diffusivity: 


86. 1 cm 2 /sec 


Horizontal Eddy Diffusivity: 


.3 
9. 2 x 10 /sec 
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Table 4. Volume and Area Data for Lake Keowee 



WATER SURFACE 

ELEVATION 
(m) (ft.) 


SURFACE AREA 
(SQ. KM.) 


STORAGE VOLUME 
CIO 9 m 3 ) 


248.9 


810 


83.3 


1.42 


245.4 


805 


78.8 


1.30 


* 243.8 


800 


74.4 


1.18 


242.3 


795 


70.1 


1.07 


240.8 


790 


65.8 


0.37 


239.3 


785 


61.2 


0.87 


237.7 


730 


56.6 


0.78 


237.1 


778 


55.1 


0.75 


** 238.2 


775 


52. 8 


0.70 


234.7 


770 


48.9 


0.62 


231.6 


760 


41.5 


0.48 


228.6 


750 


35.3 


0.36 


225.6 


740 


27.8 


0.27 


222.5 


730 


22.5 


0.19 


219.4 


720 


17.8 


0.13 


216.4 


710 


13.2 


0.08 


213.4 


700 


9.6 


0.05 


210.3 


690 


5.7 


0.02 


207.3 


680 


3.1 


0.01 


204.2 


670 


1.3 


0.004 


201.2 


660 


0.6 


0.001 


198.1 


650 


0.3 





* FuH Pond 
** Maximum Allowable Drawdown 
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Table 5. Inflows and Outflow* to Lake Keowae, August 24. 1978 



TIME 


OCONEE 
DISCHARGE 


OCONEE 
DISCHARGE 
TEMP. 


NET JOCASSEE 
FLOW 


KEOWEE 
HYDRO FLOW 


Aug. 24 
1978 


(m 3 /rain) 


°C 


(CFS) 


(CFS) 


12.00 AM 


8172. 9 


31.8 


-2638 


6638 


01.00 


8177.5 


31.8 


-6059 


6696 


02.00 


8178.6 


31.9 


-7982 


5416 


03.00 


8180.1 


32.1 


-12036 


5418 


04.00 


8185.4 


32.3 


-12066 


6756 


05.00 


8179. 4 


32.4 


-11985 


2276 


06.00 


8181.7 


32.5 


-10598 


48 


07.00 


8313.* 


32.3 


-5970 


48 


08.00 


8183.5 


32.4 


-995 


48 


09.00 


8176.4 


32.5 


100 


48 


10.00 


8183.9 


32.5 


100 


48 


11.00 


8177.1 


32.5 


100 


48 


12.00 PM 


8159.9 


32.6 


100 


48 


1.00 


8155.9 


32.6 


100 


48 


2.00 


8247.2 


32.7 


100 


48 


3.00 


8169.5 


32.7 


100 


48 


4.00 


8438.7 


32.4 


100 


48 


5.00 


7935.2 


30.7 


8726 


14224 


6.00 


75 ' 


31.0 


10905 


8432 


7.00 


8091.- 


30.3 


9823 


4544 


8.00 


8119.6 


30.2 


4557 


18 


9.00 


81 32. 8 


30.2 


6482 


48 


10.00 


8141.9 


30.2 


4183 


48 


11.00 


8143.8 


30.2 


100 


48 


12.00 AM 


8151.8 


30.3 


100 


48 
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Table 6. Inflowi and Outflows to Lake Keowee, August 25, 1978 



TIME 


OCONEE 
DISCHARGE 


OCONEE 
DISCHARGE 
TEMP. 


NET JOCASSEE 
FLOW 


KEOWEE 
HYDRO FLOW 


Aug. 25 

1978 


(ra 3 /min) 


°C 


(CFS) 


(CFS) 


12.00 AM 


8151.8 


30.3 


100 


48 


1.00 


7952. 3 


28.9 


4933 


3578 


2.00 


7097. 2 


28. 5 


838 


608 


3.00 


7098.4 


29.1 


100 


48 


4.00 


7098.8 


29.6 


100 


48 


5.00 


7103. 7 


30.0 


100 


48 


6.00 


7102.6 


30.5 


100 


48 


7.00 


7104.4 


30.7 


100 


48 


8.00 


7279.7 


30.9 


100 


3484 


9.00 


8012.8 


30.5 


100 


8624 


10.00 


7760.4 


30.7 


100 


9352 


11.00 


7759. 2 


30.9 


100 


1392 


12.00 PM 


8059.8 


30.7 


1632 


2564 


1.00 


8130.6 


30.8 


5442 


8820 


2.00 


8112.8 


30.9 


11553 


9500 


3.0C 


8141.9 


31.1 


19048 


7624 


4.0C 


8102.2 


' 31.3 


17180 


7856 


5.00 


8108.6 


31.3 


10615 


8852 


6.00 


8107.5 


31.4 


10158 


8172 


7.00 


8110.1 


31.5 


10076 


2660 


8.00 


8114.7 


31.5 


2185 


48 


9.00 


8126.8 


31.5 


29.83 


48 


10.00 


8148.3 


31.5 


3151 


48 


11.00 


8141.9 


30.4 


2647 


48 


12.00 AM 


8141.9 


30.3 

1 


100 


48 



83 



Table 7. Meteorological Data for Lake Keowee, August 24, 1978 



Time 
(Hra. from 
midnight) 


Wind Speed 
(m/s) 


Air Temp 


Dew Pt. 
Temp. 
(°C) 


Solar 

Radiation 

(w/tn») 


Wind 
Direction 


1 


1.788 


21.9 


21.94 


0.0 


30. 


2 


2.76 


21.4 


21.39 


0.0 


10 


3 


1.96 


20.8 


20.83 


0.0 


30.0 


4 


2.53 


20.8 


20.83 


0.0 


20 


5 


7.45 


20.3 


20.28 


0.0 


35.0 


6 


2.27 


20.0 


20.0 


0.0 


20.0 


7 


1.38 


19.4 


19.4 


34.85 


30.0 


8 


2.85 


20.4 


20.0 


170. 77 


15,0 


9 


1.56 


21.4 


20.83 


328. 43 


10.0 


10 


1.29 


23.3 


21.67 


592. 45 


40.0 


11 


1.69 


26.1 


23.06 


710. 94 


20.0 


12 


1.96 


27.2 


23.89 


815. 49 


30.0 


13 


2.90 


28.1 


24.44 


864. 28 


40.0 


14 


2.86 


28.9 


25.56 


838.4 


15.0 


IS 


2.90 


29.4 


25.83 


752.76 


85.0 


18 


1.82 


29.4 


25. *.8 


599. 42 


60.0 


17 


1.78 


29.4 


25.28 


397. 29 


70.0 


13 


1.42 


29.4 


26.67 


214.9 


85.0 


19 


1.47 


29.9 


26.53 


81.32 


75.0 


20 


1.16 


27.2 


25.97 


5.81 


90.0 


21 


0.71 


25.8 


25.28 


0.0 


30.0 


22 


1.82 


25.28 


25.0 


0.0 


20.0 


23 


1.16 


24.17 


24.17 


0.0 


5.0 


24 


1.81 


23.89 


23.89 


0.0 


25.0 
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Table 8. Meteorological Data for Lake Keowet, Auguet 23, 1978 



Time 
(Hrs. from 
tnidnlffht} 


Wind Speed 
(m/fj 


Air Temp 
(°C) 


Dew Pt. 
Temp. 


Solar 

Radiation 

(w/m a ) 


Hind 

Direction 

(Degreee) 


1 


1.028 


23.1 


23.0 


0.0 


55.0 


2 


2.951 


23.1 


23.0 


0.0 


40.0 


3 


2.414 


22.2 


22.2 


0.0 


20.0 


4 


2.995 


22.2 


22.2 


0.0 


55.0 


5 


1.833 


21.9 


21.9 


0.0 


30.0 


6 


1.833 


21.7 


21.7 


0.0 


35.0 


7 


2.235 


21.4 


21.4 


0.0 


25.0 


8 


1.922 


21.7 


21.7 


27.89 


40.0 


9 


1.833 


23.3 


20.0 


168.52 


10.0 


10 


0.780 


25.3 


17.2 


348. 67 


15.0 


11 


1.162 


28.9 


18.3 


534. 62 


10.0 


12 


1.565 


28.6 


16.1 


691.52 


45.0 


13 


1.654 


29.7 


16.9 


796.12 


55.0 


14 


1.967 


30.6 


16.7 


1041.35 


55.0 


15 


2.503 


30.8 


16.1 


807.74 


70.0 


18 


2.012 


30.8 


16.1 


708.96 


45.0 


17 


2.056 


30.8 


16.7 


575.30 


55.0 


18 


1.743 


30.0 


20.0 


284.74 


90.0 


19 


1.967 


29.4 


22.7 


185.96 


80.0 


20 


1.654 


28.3 


21.1 


58.11 


70.0 


21 


2.414 


27.2 


25.*" 


0.0 


65.0 


22 


1.341 


25.8 


25.6 


0.0 


90.0 


23 


1.073 


25.6 




0.0 


10.0 


24 


3.085 

i 


25.6 


25.6 


0.0 


40.0 
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Table 9. Inflows and Outflow! to Lake Keowee, February 27, 1979 



TIME 

Feb. 27, 
1978 


OCONEE 
DISCHARGE 

(m 3 /min) 


OCONEE 
DISCHARGE 
TEMP. 
C°C) 


NET JOCASSEE 
FLOW 

(C.F.S.) 


KEOWEE HYDRO 
FLOW 

(C.F.S.) 


12.00 AM 


7505. 3 


18.6 


-14395 


48 


1.00 


7498.1 


18.5 


-18754 


48 


2.00 


7492.0 


18.4 


-18805 


48 


3.00 


7492.0 


18.5 


-18713 


48 


4.00 


7491,8 


18.3 


-18698 


48 


5.00 


7494. 3 


18.3 


-18688 


48 


8.00 


7488. 2 


18.3 


-15939 


48 


7.00 


7481. 8 


18.2 


3484 


3668 


8.00 


7485. 6 


18.3 


16823 


17540 


9.00 


7488. 2 


18.2 


13503 


8488 


10.00 


7497.7 


18.3 


5470 


"096 


11.00 


7504.1 


18.3 


100 


2680 


12.00 PM 


7503.4 


18.4 


100 


48 


1.00 


7506.0 


18.5 


100 


48 


2.00 


7506.4 


18.5 


100 


48 


3.00 


7503.4 


18.5 


100 


48 


4.00 


7501.8 


18.4 


100 


48 


5.00 


7507.5 


18.4 


100 


48 


6.00 


7511.0 


18.4 


100 


48 


7.00 


7516.2 


18.4 


100 


48 


8.00 


7518.9 


18.3 


100 


'« 


3.00 


7520. 4 


18.3 


100 


48 


10.00 


7516.6 


18.2 


100 


48 


11.00 


7509.4 


18.2 


100 


48 


12.00 AM 


7507.2 


18.2 


-4382 


48 
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Table 10. Inflows and Outflows to Lake Keowee, February 28, 1979 



TIME 1 


OCONEE 
DISCHARGE 


OCONEE 
DISCHARGE 
TEMP. 


NET JOCASSEE 
FLOW 


KEOWEE HYDRO 
FLOW 


Feb. 28, 
1978 


(m 3 /min) 


(°C) 


(C.F.S.) 


(C.F.S.) 


12.00 AM 


7507. 2 


18.2 


-4382 


48 


1.00 


7616. 4 


19.8 


-18050 


48 


2.00 


7655. 5 


19.5 


-18621 


48 


3.00 


7356. 3 


19.8 


-18602 


48 


4.00 


7233. 4 


20.1 


-18692 


48 


5. GO 


7494.7 


19.5 


-18559 


48 


6.00 


7605. 3 


19.1 


-15139 


48 


7.00 


7618. 6 


19.1 


3045 


48 


8.00 


7300. 2 


19.5 


19245 


5788 


9.00 


7239. 8 


19.5 


8338 


48 


10.00 


7241. 9 


19.4 


100 


48 


11.00 


7243.0 


19.3 


100 


48 


12.00 PM 


7246. 2 


19.1 


100 


48 


1.00 


7239.0 


19.0 


100 


48 


2.00 


7241.0 


18.9 


100 


48 


3.00 


7235. 3 


18.9 


100 


48 


4.00 


7207.8 


18.9 


100 


48 


5.00 


7196.3 


18.7 


100 


48 


6.00 


7200.0 


18.9 


100 


48 


7.00 


7211.5 


19.2 


100 


48 


8.00 


7212. 


19.0 


100 


48 


9.00 


6829.4 


18.3 


100 


48 


10.00 


7189.1 


18.5 


100 


48 


11.00 


7287.8 


18.1 


-1544 


48 


12.00 AM 


7036.7 


15.8 


-15141 


48 
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Table 11. Meteorological Data for Lake Keowee, February 27. 1979 



Time 

(Hrs from 
midnight) 


Wind Speed 
(m/s) 


Air Temp 
C°C) 


1 Dew Pt. 
Temp. 
f°C) 


Solar 

Radiation 

(w/m*) 


Wind 

Direction 

(Degrees) 


1 


1.833 


-0.33 


-2.78 


0.0 


15° 


2 


1.073 


-0.72 


-1.67 


0.0 


75° 


3 


2. 325 


-1.61 


-1.61 


0.0 


60° 


4 


. 1. 565 


-2.22 


-2.28 


0.0 


15° 


5 


2.056 


-1.83 


-1.89 


0.0 


50° 


6 


1.788 


-2.17 


-2.22 


0.0 


85° 


7 


2.012 


-2.72 


-2.78 


20.94 


85° 


8 


2.280 


-1.67 


-2.78 


195.39 


60° 


9 


0.626 


0.01 


-3.33 


369.85 


5° 


10 


1.386 


3.06 


-2.22 


544. 31 


75° 


11 


1.609 


5.83 


-2.22 


855. 96 


15° 


12 


1.788 


8.83 


-1.39 


725. 75 


40° 


13 


3.129 


11.06 


-2.78 


746.68 


80° 


14 


2.593 


12.28 


-5.0 


704.81 


70° 


15 


1.520 


13.39 


-5.56 


579.20 


80° 


16 


1.207 


13.89 


-5.56 


383.81 


75° 


17 


1.565 


13.83 


-5.61 


146.55 


55° 


18 


1.609 


13.72 


-3.33 


20.94 


15° 


19 


2.056 


11.72 


-4.44 


0.0 


30° 


20 


1.162 


S.72 


-2.78 


0.0 


25° 


21 


2.772 


8.33 


5.28 


0.0 


55° 


22 


2.861 


7.78 


5.56 


0.0 


55° 


23 


2.995 


7.00 


5.28 


0.0 


50° 


24 


1.386 


5.28 J 


3.89 


0.0 


60° 
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• Table 


12. Meteorological Data for Lake Keowee, February 28, 1979 


Time 

(Hrs. iron 
mid night) 


Wind Speed 
Cm/s) 


Air Temp 
(°C) 


Dew Pt. 

Temp. 

(°C) 


Solar 
Radiation 

(w/m a ) 


Wind 

Direction 

(Degrees) 


1 


2.280 


4.94 


4.17 


0.0 


55 


2 


2.101 


4.00 


4.17 


0.0 


55 


3 


1.967 


3.00 


2.22 


0.0 


60 


4 


2.414 


2.28 


1.67 


0.0 


50 


5 


2.593 


2.61 


2.22 


0.0 


80 


6 


4.068 


4.11 


3.06 


0.0 


75 


7 


3.398 


3.83 


0.28 


20.94 


60 


8 


3.934 


4.22 


-2.78 


153.52 


60 


9 


4.515 


6.67 


-3.89 


279.13 


85 


10 


4.381 


7.78 


-5.00 


439.64 


60 


11 


3.710 


9.28 


-5.00 


676.90 


70 


12 


2.235 


10.78 


-3.89 


907.18 


45 


13 


1.878 


11.72 


-3.89 


537. 33 


60 


14 


1.609 


12.28 


-3.89 


565. 24 


50 


15 


3.442 


13.33 


-3.33 


690.85 


45 


16 


2.950 


13.61 


-2.78 


0.0 


45 


17 


2.325 


13.33 


-1.67 


0.0 


45 


18 


2.012 


12.72 


-0.56 


0.0 


5 


19 


2.235 


11.39 


1.39 


0.0 


65 


20 


2.950 


10.39 


2.28 


0.0 


55 


21 


2.772 


10.61 


2.22 


0.0 


85 


22 


2.638 


10.44 


2.22 


0.0 


50 


23 


2.191 


10.00 


2.22 


0.0 


65 


24 


2.369 


9.83 


2.22 


0.0 


70 
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Table 13. Catalogue of Runs for 3-D Lake Keowee 



S 



Run 
Identification 


Wind 
(M/Sec) 


Current 
(Cm/Sec) 


Discharge 
Velocity 
(cm/Sec) 


Discharge 
Temperature 


Depth (Meters) 


Remarks 


Variable 


Constant 


L001 
(Selection of 
B.C. '» and I.C.'s 


Mo 


No 


5.65 


32.3 


NO 


Yes (16 
Meters 


Total Simulated 
Timet 21.60 Hrs. 


LOO 2 
(Model Execution 
for One physical 
parameter - Vari- 
able Depth) 


No 


No 


5.65 


32.3 


Yes 


NO 


Total Simulated 
Timet 8.64 Hrs. 


L003 
(Model Execution 
for Two Physical 
parameters - 
Depth and Wind 


4.00 


No 


5.65 


32.3 


Yes 


NO 


Total Simulated 
Timet 8.64 Hrs. 


L004 
(First Execution 
for Archival Data 


12.83 


0.69 


5.65 


32.4 


Yes 


NO 


Total Simulated 
Timet 8.64 Hrs. 


LOOS 
(Second Execution 
for Archival Data) 
Up Date: L006 
(For Ground Truth 
Data) 


12.83 


0.69 


5.65 


32.4 


No 


Yes (16 
Meters) 


Total Simulated 
Timet 32.4 Hrs. 


L006 
August 24/25 Data 
Base Comparison 


variable 
(0.76 to 
3.09) 


variable 
(0.61 
to 1.12) 


7.42 


18.4 


No 


Yes (16 
Meters) 


Total Simulated 
Timet 48 Hours 


LOO 7 
February 27/28 
Data Base Com- 
parison 


variable 

(1.61 
to 4.52] 


Variable 

(0.61 
to 4.14) 


6.84 


18.4 

— 


NO 


Yes (16 
1 Meters) 


Total Simulated 
Timet 48 Hours 



Table 14. Root Mean Square Difference Between IR 
and Predicted Temperatures 



Time 


RMS Difference 


Morning, August 24, 1978 


0.55 °C 


Morning, August 25, 1973 


0.34 °C 


Afternoon. February 27, 1979 


0.82 °C 


Morning, February 28, 1979 


0.01 °C 
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